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METHOD FOR REDUCING SKEW 
IN A REAL-TIME CENTROID CALCULATION 

Field of the Invention 

The present invention relates to real-time signal 
processing, and more particularly to the real-time determination 
of a center value (centroid) of a waveform. 

Background of the Invention 

Often it is necessary to determine precisely a center value 
of a waveform, such as the wavelength value at which a peak of a 
waveform occurs in a signal representing light energy over a 
range of wavelengths, and to do so in a way that involves a 
small computational burden, so as to perform the calculation in 
"real-time" (i.e. performing the calculation at a rate 
sufficient to keep pace with the arrival of each new set of 
calculation inputs) . For example, in many situations, it is 
necessary to determined the center value of a waveform in a 
signal of light reflected from a Bragg grating inscribed in an 
optical fiber. An optical fiber in which a Bragg grating is 
inscribed can be used as a component of a pressure sensor or as 
a component of an optical add/drop multiplexer. The wavelength 
of the light reflected by the Bragg grating (i.e. the center 
value of the peak in a signal corresponding to the reflected 
light) conveys information of use in the sensor and add/drop 
multiplexer applications. 

The precision required for different applications varies, 
but a precision of better than 1 part in 300 of the full width 
at half maximum (FWHM) of a waveform is often required. A 
simple technique to determine the center of a waveform (i.e. 
e.g. the position of a peak in the waveform) uses a centroid (or 
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center of mass) calculation over an area of interest. For an 
evenly sampled waveform having n sampled points, the basic 
equation for the calculation of the centroid of the waveform is: 

K=- n (i) 

where Ai is the sampled amplitude at wavelength Vi (or some other 
parameter, besides wavelength) . As seen in equation (1) , each 
sample point is weighted by the amplitude at the sample point. 

A centroid calculation according to equation (1) can be 
performed quickly is therefore useful as a real-time computer 
algorithm. Equation (1) can be used for example to determine 
the center wavelength V c of a waveform where n measurements Ai 
of power level are taken at n different wavelengths Vi. 

The number n of sampled points in the calculation is 
critical when using equation (1) . Ideally, the samples span the 
entire waveform of interest. However, in practical applications 
this is not always possible due to factors such as sampling 
rates and memory requirements. When the samples do not span the 
entire waveform, a fixed number of sample points around the peak 
value of a waveform are used. However, the locations of the 
sampled points affect the determination of the waveform center, 
especially when the points are not perfectly symmetric about the 
actual waveform peak. 

In a typical real-time application in which a centroid is 
determined, only a small number of sampled points of a waveform 
are available from which to determine the centroid of the 
waveform, sometimes as few as seven points and sometimes even 
fewer. Even with such a small number of sampled points however, 
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if the sampled points are symmetric with respect to the true 
centroid, the calculation according to equation (1) will provide 
the true centroid. On the other hand, if the sampled points are 
unevenly distributed with respect to the true centroid, the 
5 calculation according to equation (1) will in general be skewed. 

(The use of equation (1) for some uneven distributions could 
still yield the correct centroid, as long as the sum of the 
moment arms (amplitude times wavelength) on one side of the true 
centroid balances the corresponding sum on the other side.) 

10 A co-pending, co-owned application serial number 087447, 



filed May 29, 1998, entitled METHOD FOR IMPROVING THE ACCURACY 
IN THE DETERMINATION OF A WAVEFORM CENTER, discloses a method of 
correcting for skew in a calculation of a centroid of a waveform 
based on eliminating some of the sample points used in the 
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calculation according to equation (1) . In some situations, 
especially when the number of sample points is already a very 
small number, ignoring some of the sample points tends to 
increase error in the calculation of the centroid. 
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What is needed is a method for reducing skew in the 
cal culation of the centroid of a waveform not based on only 
eliminating any of the sampled points, but rather based on 
adding points to be used in the calculation of the centroid. 



Summary Of The Invention 
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Accordingly, the present invention provides a method and 
corresponding apparatus for determining the centroid (V c ) of a 
waveform signal indicating how the output of a system, such as 
the light reflected by a Bragg grating inscribed in an optical 
fiber, responds to an input signal, such as incident light, as 
the value of a characteristic of the input signal, such as the 
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wavelength of the incident light, is varied over a predetermined 
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range, the waveform being sampled at a set of parameter values 
(Vi, i=l, . n) yielding a corresponding set of sampled 
amplitudes (At, i = l, . . . , n) 7 each parameter value and 
corresponding amplitude forming a sampled point (Vi, Aj.) , the 
method including the steps of: selecting an amplitude at which 
to create an interpolated point; interpolating a first parameter 
value corresponding to the amplitude selected in the step of 
selecting an amplitude; and performing a centroid calculation 
using only the sampled points with an amplitude greater than a 
predetermined threshold. 

In a further aspect of the invention, the amplitude 
selected in the step of selecting an amplitude is less than 
approximately twenty per cent of the maximum sampled amplitude. 

In another further aspect of the invention, the centroid 
(V c ) is calculated using as a formula: 



In yet another further aspect of the invention, the 
waveform is sampled in the presence of background noise, and the 
method also includes the steps of: estimating the background 
(Bi) for each value in the set of parameter values at which 
sampling is performed; and reducing the amplitude (Aj.) of each 
sampled amplitude by the background (Bi) so estimated. 

In still yet another further aspect of the invention, there 
is, among the sampled amplitudes, a maximum sampled amplitude, 
and, in addition, the method also includes the step of 
interpolating a second parameter value to correspond to the 
amplitude selected in the step of selecting an amplitude, the 
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second value on the opposite side from the first interpolated 
value of the maximum sampled amplitude. 

Brief Description of the Drawings 

The above and other objects, features and advantages of the 
invention will become apparent from a consideration of the 
subsequent detailed description presented in connection with 
accompanying drawings, in which: 

Fig. 1 is a graph of the sampled points of a waveform shown 
superimposed on a waveform that is sampled would yield the 
sampled points, and showing the addition of a new point to 
the set of sampled points; 

Fig. 2 is a graph of some of the sampled points of Fig. 1 
and the added sampled point, with an amplitude of each 
sampled point reduced by an amount sufficient to eliminate 
the smallest amplitude sampled point of the original set of 
sampled points; 

Fig. 3 is a graph of all of the original sampled points of 
Fig. 1, and showing the addition of two new points to be used 
in the calculation of the centroid of the waveform; 

Fig 4 is a graph of some of the sampled points of Fig. 3, 
as well as the points added to the original set of sampled 
points, with the amplitude of each of the sampled points and 
that of the added points reduced by an amount sufficient to 
eliminate use of several of the original sampled points in 
the calculation of the centroid of the waveform; and 

Fig. 5 is a flow chart of the method of the present 
invention. 
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Best Mode For Carrying Out The Invention 

In order to reduce skew in a real-time calculation of the 
centroid V c of a waveform sampled at a set of wavelengths Vi (or 
some other physical parameters on which the waveform depends) 
and so observed to have corresponding amplitudes Ai, the present 
invention involves adding to the set of sample points 
(wavelength-amplitude pairs) of the waveform either one or two 
interpolated points having an amplitude equal to a predetermined 
amplitude that is typically a small fraction (usually 2 0% or 
less) of the maximum sampled amplitude. Whether one or two 
points are added depends on whether an existing sample point has 
an amplitude approximately equal to the predetermined amplitude. 

Referring now to Fig. 5, the method begins with a step 51 
of providing a predetermined amplitude at which one or two 
points are to be interpolated. Next, in a step 52, a search of 
the existing sample points is made to determine whether any such 
point has an amplitude approximately equal to the predetermined 
amplitude and if so to identify such a point. 

Adding only a single interpolated point 

Referring still to Fig. 5, if the sampled points include a 
point having an amplitude that is approximately equal to the 
predetermined amplitude, then only a single interpolated point 
is added to the set of points used in the centroid calculation. 
But first, in a step 54 (optional), the predetermined amplitude 
is adjusted to equal more precisely the amplitude of the 
identified existing sample point. 

Referring now to Figs. 1, 2 and still to Fig. 5, the method 
of the present invention in case of adding only a single 
interpolated point is shown as including the steps of: (in a 
step 55) creating the interpolated point Nl, at a predetermined 
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amplitude, to be used in combination with some of the sampled 
points PI, P2, . .., P7; (in a step 56) subtracting from each 
amplitude an offset (the same for each amplitude) sufficient to 
cause the interpolated point to have a small non-zero amplitude, 
but also sufficient to cause the sampled points outside the 
interpolated point (PI in Fig. 1) to have a negative value (see 
Fig. 2) ; and (in a step 57) performing the centroid calculation 
using one or another centroid formula, such as that given by 
equat ion ( 1 ) , i.e. 

' c n ' 

but including in the calculation only the originally sampled 
points having a positive amplitude after subtracting the 
predetermined constant offset, (i.e. points P2 , P7 of Fig. 

1) as well as the interpolated point Nl, as indicated in Fig. 2. 

Another example of a formula that can be used to calculate 
the centroid, instead of that given by equation (1) , is given by 
equation (5) below. 

In a sense what is being done here is a reverse 
interpolation, in that the ordinate (amplitude) or y-value of 
the interpolated point is given (predetermined so as to have the 
interpolated point become an endpoint) and its abscissa 
(wavelength) or x-value is to be determined based on some method 
of interpolation, such as linear interpolation or cubic spline 
interpolation, or based on curve fitting using an assumed form 
of the equation for the waveform is known. 

Still referring to Figs. 1, 2 and 5, the interpolated point 
Nl is intended to even out (relative to the centroid) the 
distribution of wavelength-amplitude pairs used in the 
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calculation of the centroid. Ideally, the resulting 
distribution is symmetric about the true centroid. In such a 
case, the calculation according to equation (1) will yield the 
actual centroid value, i.e. the wavelength at which the actual 
maximum of the waveform occurs . 

Still referring to Figs. 1, 2 and 5, according to the 
invention, the interpolated point that is added is to become an 
endpoint of the set of points used in the calculation of the 
centroid, as indicated in Fig. 2, and is therefore selected (in 
a step 51) to have a predetermined amplitude that is typically a 
small fraction (usually less than 20%) of the amplitude of the 
sampled point having the maximum amplitude (i.e. the amplitude 
of P4 in the case of the waveform of Fig. 1) . 

There is computational burden associated with adding an 
interpolated point. In some cases, as illustrated in Fig. 1, 
there is an existing sample point (P7) having an amplitude at or 
near the (predetermined) amplitude at which an interpolated 
point (Nl) is to be added. In such a case, the predetermined 
amplitude is adjusted to equal the amplitude of the existing 
sampled point, and then only a single interpolation need be 
performed to add one point on the opposite side of the peak from 
the existing sampled point, and at the same amplitude as the 
existing sampled point. Then the offset is determined and 
subtracted from each sampled amplitude as well as from the 
amplitude of the interpolated point so that the interpolated 
point (Nl) and the existing sample point (P7) become the 
endpoints of the waveform, with each having the same small, non- 
zero amplitude. Alternatively, the offset can be applied before 
the interpolation. 

There is a significant saving of computational burden in 
only having to add a single interpolated point. The burden of 
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adding one interpolated point, as illustrated in Fig. 1 where an 
existing sampled point is found to have approximately the 
predetermined amplitude, is approximately only one half that for 
adding two interpolated points, because it turns out to be 
5 possible to determine whether an existing sampled point has an 

amplitude approximately equal to the predetermined amplitude in 
as little as a single clock cycle. 

Referring now in particular to Fig. 2, the same waveform is 
shown in Fig. 2 as in Fig. 1 except reduced in amplitude by the 
10 predetermined offset. The result is that point PI has a 



negative amplitude and is therefore discarded (not used in the 
calculation of the cent roid) . The points used in the 
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calculation of the centroid after subtracting the predetermined 
offset then include three points on either side of the maximum 
sampled amplitude P4 . 
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Ideally an estimated background is also subtracted from 
each sampled amplitude, before subtracting the constant offset, 
the estimated background being in general different for 
different wavelength values. A subtraction of a constant offset 
from all sample points and also from the interpolated point has 
no effect on the calculation of the centroid, except that it 
drives some of the sampled amplitudes to a negative value 
(specifically the sampled amplitude for point PI in Fig. 1) , so 
that such points are not then used in the calculation of the 
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centroid. On the other hand, the subtraction of a variable 
background can reduce skew. 



Creating two interpolated points 
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Referring now to Fig. 3 and Fig. 5, the method of the 
present invention is shown in case there is no existing sample 
point having an amplitude approximately equal to the 
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predetermined amplitude so that two interpolated points are 
created. Fig. 3 shows the same waveform as in Fig. 1, with the 
same sampled points PI P2 ... P7, but indicates the addition in 
a step 53 of two interpolated points Nl N2 on either side of the 
maximum sampled amplitude P4 . As explained above, the addition 
of two interpolated points involves nearly twice the 
computational burden as the addition of one interpolated point 
at an amplitude for which there is already a sampled point. 

Referring now to Fig. 4, the waveform of Fig. 3 with the 
sampled points P3 ... P6 inside the two interpolated points Nl 
N2 shows the reduction of amplitude (step 56) at each of the 
sampled points at an amount sufficient to eliminate all sample 
points outside of the two added interpolated points Nl N2 . 

Scope of the invention 

It is to be understood that the above-described 
arrangements are only illustrative of the application of the 
principles of the present invention. In particular, it is 
obvious that it is not actually necessary to subtract a 
predetermined offset from the amplitude of each of the sampled 
and interpolated points. It is enough to determine whether a 
given point has an amplitude greater than the offset and to 
include it in the calculation only in such a case, so that the 
offset serves as a predetermined threshold. Numerous other 
modifications and alternative arrangements may be devised by 
those skilled in the art without departing from the spirit and 
scope of the present invention, and the appended claims are 
intended to cover such modifications and arrangements. 
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